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The Poisson-Dirichlet distribution arises in many different areas. 
The parameter 9 in the distribution is the scaled mutation rate of a 
population in the context of population genetics. The limiting case of 
6 approaching infinity is practically motivated and has led to new, in- 
teresting mathematical structures. Laws of large numbers, fluctuation 
theorems and large-deviation results have been established. In this 
paper, moderate-deviation principles are established for the Poisson- 
Dirichlet distribution, the GEM distribution, the homozygosity, and 
the Dirichlet process when the parameter 9 approaches infinity. These 
results, combined with earlier work, not only provide a relatively 
complete picture of the asymptotic behavior of the Poisson-Dirichlet 
distribution for large 6, but also lead to a better understanding of 
the large deviation problem associated with the scaled homozygos- 
ity. They also reveal some new structures that are not observed in 
existing large-deviation results. 

1. Introduction. For ^ > 0, let cri(^) > o"2(^) > ••• be the points of a 
nonhomogeneous Poisson process with mean measure density 

and (T{9) = J2°Zi(^i{d)- Set 



(1.1) P{9) = iP,{9),P2{e), 



Then it is known that P(^) and cr{9) are independent, and cr{9) is a Gamma(0, 1) 
distributed random variable. The law of P(^) is called the Poisson-Dirichlet 
distribution with parameter 9, and is denoted by PD(0). 
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Let Uk, k = 1,2, . . . , be a sequence of i.i.d. random variables with common 
distribution, Beta(l,^). Set 

(1.2) Xi{e)=Ul, Xnie) = il-Ul)---{l-Un^l)Un, n>2. 

Then with probabihty one 

oo 
k=l 

and the law of {Xi{9),X2{9), . . .) is called the GEM distribution, denoted by 
GEM(6l). The law of the descending order statistics X^i){9) > X(2)(6') > • • • 
of Xi{e),X2i9),... is also FD{9). 

Let ^fc, = 1, . . . , be a sequence of i.i.d. random variables, independent of 
P(^), with a common diffusive distribution on [0,1], that is, ^{{x}) = 
for every x in [0, 1]. Set 

oo 

(1-3) Eg,, = Y,Pk{9)6^^. 

k=l 

We call the law of '^e^u, the Dirichlet process, denoted by Dirichlet(0, z^). 

The Poisson-Dirichlet distribution was introduced by Kingman [21] to 
describe the distribution of gene frequencies in a large neutral population 
at a particular locus. The component Pk{9) represents the proportion of 
the kth. most frequent allele. If e is the individual mutation rate and N^. is 
the effective population size, then the parameter 9 = 4:Nee is the scaled 
population mutation rate. The GEM distribution can be obtained from 
the Poisson-Dirichlet distribution through a procedure called size-biased 
sampling. It provides an effective way of doing calculations involving the 
Poisson-Dirichlet distribution. The name, GEM distribution, was coined by 
Ewens after Grifffiths, Engen and McCloskey for their contributions to the 
development of the structure. The Dirichlet process first appeared in [11] 
in the context of Bayesian statistics. It can be viewed as a labelled version 
of the Poisson-Dirichlet distribution. More background information can be 
found in [8]. 

For any integer m > 2, consider a random sample of size m from a pop- 
ulation following the Poisson-Dirichlet distribution. Given the population 
proportion, p = {pi,P2, ■ ■ ■), the probability that all samples are of the same 
type is given by 

oo 

i=l 

The quantity i^2(p) is called the population homozygosity. It is an im- 
portant statistic in population genetics. For general ni, we refer to Hm.{p), 
as the homozygosity of order m. 
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Consider a family of random variables {1a : A > 0}. Assume a law of large 
numbers holds; that is, Yx converges in distribution to a constant c as A 
approaches infinity. A fluctuation theorem such as the central limit theorem 
is a statement that there exists a function 5(A) approaching infinity for large 
A such that 

b{X){Yx-c)^Y, A^oo, 

where y is a nontrivial random variable and "=^>" denotes convergence in 
distribution. A large-deviation result is concerned with estimates of proba- 
bilities P{Yx — c S ^4} for measurable sets A. A moderate-deviation result 
lies between the fluctuation theorems and large deviations. It is concerned 
with estimates of probabilities P{a{X){Yx — c) E A} for measurable sets A, 
where a(A) is an intermediate scale between 1 and b{X). 

The objective of this paper is to establish moderate-deviation principles 
(henceforth MDP) for GEM(6l), PD(6'), the homozygosity and Dirichlet(6l, v), 
when 6 approaches infinity. 

The study of the behavior of P(^) = {Pi{9), P2{0), . . .) for large 6, goes 
back to the seventies. In Watterson and Guess [29], E[Pi{6)] was shown to 
be asymptotically log9/9. Griffiths [15] obtained the explicit weak limit of 
9P{d) and a central limit theorem for the population homozygosity. The 
limiting case of large 9 is equivalent to a situation where the mutation rate 
per individual is fixed and the effective population size is large. Motivated 
by the work of Gillespie [12] on the role of population size in molecular 
evolution, there have been renewed interests in the asymptotic behavior of 
FB{9) for large 9 (see [4, 9, 18, 19, 20]). In particular, in [19], central limit 
theorems are obtained for the homozygosity of order m. Large deviations are 
established in [4] for PD(0) and the homozygosity, and in [9] for the GEM 
distribution. Large deviations for Dirichlet(0, z^) can be found in [3, 23]. 

Although MDP is a natural mathematical object that warrants rigorous 
investigation and our study does reveal some new mathematical phenomena, 
the real motivation for this work comes from the results in [4, 19]. It was 
shown in [19] that, as 9 goes to infinity, 

(1-4) —-H^iP{9))^l 

1 [m) 

and 

(flm— 1 \ 
^HUP{9))-lj^Z{m), 

where Z{m) is a normal random variable with mean zero and variance 

— . These are the law of large numbers and central limit theo- 
ry (m) ° 

rem for p™ s Hm{P{9)). A natural companion to these limit theorems is the 
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large deviations p^^^ Hm(P{0)) from one, or equivalently the large devia- 
tions of Yj^Hm(P{0)) — 1 from zero. Unfortunately this problem is still 
open. The large deviation principle established in [4] is concerned with the 
deviations of Hm{P{0)) from zero. The scale difference between Hm{P{9)) 
and '^Hm{V{e)) is of order of e^-\ Multiplying ^^H^{V{e)) - 1 by 
a factor O'^ , places us in the territory of MDP. One would hope that the 
study of MDP will shed light on resolving the large deviation problem which 
corresponds to 7 = 0. The MDPs we obtain require that 7 is bigger than 
a strictly positive number. Thus a gap exists between the MDPs and the 
LDP. This seems to indicate that a large deviation principle may not exist 
for ^^HminO)) - 1. 

This paper is organized as follows. The basic terminology of LDP, MDP 
and a comparison lemma are given in Section 2. In Section 3, we discuss the 
fluctuation theorems associated with PD(0), Dirichlet(0, v) and the homozy- 
gosity. A new proof is given for the central limit theorem of homozygosity 
in [19], using Campbell's theorem. A MDP for GEM is established in Sec- 
tion 4. Section 5 deals with the MDP for PD(0). Since the condition of the 
Gartner-Ellis theorem is not satisfied, we prove the result by direct calcu- 
lation. The MDP obtained in Section 6, is for the homozygosity, for which 
the MDP holds in a narrower range of scales. The proof is based mainly on 
Campbell's theorem. In the MDP literature, general results such as those in 
[7, 13, 30], usually require the finiteness of exponential moments in a small 
neighborhood of zero so that the Laplace method can be used. Here the ex- 
ponential moment is infinite on the positive half-line. One way to deal with 
the infinite exponential moment is to verify Ledoux's condition in [22]. Since 
this does not seem easy to do, we choose the truncation method instead. Fi- 
nally in Section 7, we establish the MDP for Dirichlet(0, z^). Compared to 
the Sanov theorem, the LDP rate function for Dirichlet(^, z^) is a reversed 
form of relative entropy. Here the MDP rate function for Dirichlet process 
is the same as the MDP rate function for the empirical process of an i.i.d. 
random sequence with common distribution v. When v is supported on a 
finite number of points, one can see this clearly from the fact that both the 
relative entropy and its reversed form have the same second-order derivative 
at u. 

The MDPs for the Poisson-Dirichlet distribution and GEM have a differ- 
ent speed from the MDPs for the homozygosity and the Dirichlet process, 
the latter having a more standard structure. One explanation for this is 
that in the cases of the Poisson-Dirichlet distribution and GEM, we are 
concerned with partial information such as alleles with a certain propor- 
tion size or age order, while for the homozygosity and the Dirichlet process, 
all alleles contribute. One expects that similar results and structures exist 
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for the two-parameter Poisson-Dirichlet distribution and Dirichlet process 
[10, 26]. 

2. Preliminaries. In this section we introduce the terminology on LDP 
and MDP used in this paper, and prove a comparison lemma that plays an 
important role in proving the main results. Comprehensive coverage on LDP 
techniques can be found in [6]. 

Definition 2.1. Let E hea Pohsh space with metric d, and {Yq : > 0} 
be a family of ii^- valued random variables. Denote the law of Yg hy Pq. 

(1) The family of probability measures {Pg :6 > 0} (or the family {Yg : 9 > 
0}) is said to satisfy a LDP with speed \{9) and rate function /(•), if for 
any closed set F and open set G in E 

limsupA(6l)logPe{F} < - inf I{x), 

liminf A(6')logPe{G} > - inf I{x), 
for any c> 0,{x: I{x) < c} is compact. 

In short form, we say {Pg , I {■) , X{9)) satisfies a LDP. 

(2) The family {Pg ■.9 > 0} is said to satisfy a local LDP with speed X{9) 
and rate function /(•), if for every x in E 

limlimsupX{9)logP{d{Yg,x)<5} 

= limlimini X{9) log P{d{Yg,x) < 6} = -/(x), 

and for any c> 0, {x: I{x) < c} is compact. 

(3) The family {Pg ■.9>0} is exponentially tight with speed X{9) if for 
every L > 0, there is a compact set Kl in E such that 

limsup A(^) logP{Yg i Kl} < -L. 

Remark 2.1. It is known that a local LDP combined with exponential 
tightness implies the LDP (cf. [27]). 

Definition 2.2. We use to denote convergence in distribution. 

(1) The family {Yg:9>0} is said to satisfy a fluctuation theorem if there 
exist functions b{9), c{6) and a flnite nondeterministic random variable Z 
such that 

lim b{9) = oo, b{9) [Yg - c{9)] ^ Z, 9 ^ oo. 
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(2) Assume that the family {Yg : > 0} satisfies the fluctuation theorem 
above. Let a{9) satisfy 

(2.1) lima(^) = oo, lim ^ = 0. 

The family {Pg : > 0} or equivalently the family {Yg : > 0} is said to 
satisfy a MDP with speed \{9) [depending on a{6)] and rate function /(•) 
if the family {a{9)\Yg — c{9)] :0 > 0} satisfies a LDP with speed A(^) and 
rate function /(•). Thus the MDP for {Yg:9> 0} is the LDP for {a{9)[Yg - 
c{9)]:9>{)}. 



The next lemma is a useful tool in deriving the MDPs of this paper. 



Lemma 2.1. Let {^e : > 0} and {rjg > : > 0} be two families of real- 
valued random variables. Assume that for any 5 > 

(2.2) limsupA(6l)logP(|r/e - 1| >6) = -oo. 

Then {P{^g G ■), I{x), X{9)) satisfies a LDP iff {P{CgVe G ■) , I (x) , X{9)) sat- 
isfies a LDP. 

Proof. For any 6 > 0, choose 5 = min{|, ^}. Then it is clear that 

{\Vg'-M>S}c{\ve-l\>S}, 
which, combined with (2.2), implies 

(2.3) limsupA((9)logP(|r?g ^ - 1| >5) = -oo. 

9-»oo 

For any j; e R, 7 > and 6>0, 

Pil^em - x| < 7) 

< P{\Vg' - 1| > <5) + P{\Vg' -1\<S, -x\< ^Vg' + \x\\Vg' " 1|) 

< PilVg^ - 1| > <5) + P{\Cg -x\< 7(1 +6) + \x\6) 
which implies 

lim (^^^^^ Xi9)log P{\CgVe - x| < 7) < lim (^^^^^ X{9)logP{\Cg - x| < 7). 
Symmetrically, 

lim (^'^I'P) Xi9)logP{\Cg - x| < 7) < lim (^^^^^^ Xi9) log Pd^gVe - x\ < 7). 
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Furthermore, for any L > 0, 

P{\^eve\ >L)< P{\rje - 1| > <5) + Pm > (1 + Sy^L), 

pm >L)< p{\ij^' -i\>6)+ Pi\^em\ > (1 + sy'L). 

Thus the exponential tightness of {S^e : 6* > 0} is equivalent to the expo- 
nential tightness of {S,er]g : > 0}. The lemma now follows from Remark 2.1. 
□ 



3. Fluctuation theorems. We start this section with a discussion of the 
asymptotic behavior of the random variable cr(6) for large 0. It plays a key 
role in connecting the Poisson process to the Poisson-Dirichlet distribution. 
To put our MDP results into perspective, we present in this section sev- 
eral known fluctuation theorems for the Poisson-Dirichlet distribution, the 
Dirichlet process and the homozygosity of order m. 

3.1. Gamma distribution. Recall that cr(^) is a Gamma(0, 1) random 
variable with density 



(3.1) J-u^-V", 0<n<oo, 
and exponential moment 

(3.2) i?[eMe)] = I 

[ oo, else. 

Let 

(3.3) A(t) = <(l°gT37' 

oo, else. 

Routine calculations and Gartner-Ellis theorem lead to the following the- 
orem: 



Theorem 3.1. When 9 approaches infinity, the following hold: 

(a) lime^oo ^ = 1- 

(b) Ve{^ -l)^Z, where Z is a standard normal random variable. 

(c) The family of the laws of a{9)/9 satisfies a LDP with speed 1/9 and 
rate function 



(3.4) I{u) 



u—1 — logu, ti > 
oo, else. 
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Let a{6) be a positive function satisfying 
(3.5) lim ^ = 0, lim a{6) = oo. 



Corollary 3.1. For each 5 > 0, 



(3.6) 

and 

(3.7) 



aid) ^ „ 
lim sup log P 



lim sup ^ ^ log P 



a{9) 



a{9) 



>6 



>6 



-oo 



-oo. 



Proof. Equality (3.6) is derived directly from Theorem 3.1. Since 



one gets (3.7) from (3.6). □ 



(3.8) 



Let a{6) be a positive function satisfying 

a{9) 



lim ^^^^ = 0, lim a(9) = oo. 
e-,oo ^JQ e^oo 



The following theorem is standard. 

Theorem 3.2. The family of the laws of a{9)/9 satisfies a MDP with 
speed a'^{9)/9 and rate function 

„,2 



(3.9) 



S{u) 



-oo < u < oo. 



3.2. Fluctuations. Consider a nonhomogeneous Poisson process with mean 
measure 



e "du, 



-oo <U < +00. 



Let Ci ^ C2 ^ ■ ■ ■ be the sequence of the points of the nonhomogeneous 
Poisson process in descending order. Then for each r > 1 the joint density 
of iCi, ■■■,Ct) is 

(3.10) e~^i^i'^'e~''~"\ -00 <Ur <■ ■■ <ui <oo. 

Let P{9) = log9 — log log 0. The following result is obtained in [15]. 
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Theorem 3.3. The sequence {9Pi{e) - l3{e),9P2{e) - p{9), .. .) converges 
to {Ci,C2,---) in distribution as 9 tends to infinity. 

The next theorem is obtained in [19]. We give a different proof here using 
Campbell's theorem. 

Theorem 3.4. Let 

Ak{9) = V9{^^Hk{V{9)) - l) , = 2,3, . . . 

and Ae = {A2{9),A^{9),...). Then 

(3.11) Ae^A, 6*^00, 

where A is a -valued random element and for each r>2, (A2, . . . , A^,) 
has a multivariate normal distribution with zero mean and covariance matrix 

r.i9^ r (A A\ r(fc + o-r(fc + i)r(/ + i) 

(3.12) Cov{Ak,Ai) = r{k)r{l) ' k,l = 2,...,r. 

Proof. For each A; > 1, set 

B, = {Bi{6l),...). 

For each fixed r > 1 and any (ai, . . . , a,.) in , set 

It follows from Campbell's theorem that 

_E(e[**ELi°fe'^fcW]) 

(3.13) =e[-^*SI=i"^^lexp|0^°°(e**^(^) -l)y-ie-^dy| 



exp< 



t^ ^ ru + k) 



Let B = (Bi, . . .) be such that for each r > 1, {Bi, . . . ,Br) is a multivariate 
normal random vector with zero mean and covariance matrix 

(3.14) I|±|, ,.. = 1,...,. 
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Then (3.13) implies that Bg converges in distribution to B. 

For r > 2, it follows from (1.1) that the following relation holds between 
{A2{e),...,Ar{e)) and (^2(0), . . . , i3.(^)): 




It follows from the convergence of Bg to B that 



(3.16) 




in distribution. 



By Theorem 3.1 and basic algebra, one gets 

(3.17) v5((_|_)''_i)^_,B, 

By (3.15)-(3.17), one gets 

r r 

(3.18) Y.'^kAk{e)^Y.'^k{Bk-kBi). 

k=2 k=2 

The theorem now follows from the fact that the covariance of {B}. — kBi) 

and(i?,-/i?Ois^i^^±^^Tlffi?^- ° 

Let {X[t),t G [0,00)} be a Gamma process; that is, a stochastic process 
with stationary independent increments and right-continuous paths with 
X(0) = and such that X[l) has an exponential distribution with parameter 
1. For each Borel measurable set A, define 

XeAA) = X{eu{A)) 

and 

x{ev{A)) 

Ze,M)- x^Q^ , 

where ;y is a diffusive distribution on [0, 1], that is, v{{x}) = for every x in 
[0,1]. Set i/(t) = z^([0,t]), Xe,,{t)=Xe,,{[Q,t\) and Ze^,{t) = ^^,,([0,*]). Then 
Ze.ui-)-, as a random measure, is distributed as DirichIet(0, z^). Let D([0,1]) 
be the space of all real- valued cadlag functions defined on [0, 1] that are left 
continuous at 1, equipped with the topology of uniform convergence. Then 
the functional central limit theorem for processes with independent incre- 
ments yields immediately that {Xg ,j{t) — 6u{t)) / \/6 converges to B{u{t)) in 
distribution on D([0, 1]), where B(t) is a standard Brownian motion (cf. VII 
3.5 in [16], page 373). This, combined with the fact that X{6)/9 converges 
to 1, implies the following result: 

Theorem 3.5. The family of processes {\^{Zg^^{t) — v{t)),9 > 0} con- 
verges to B[u{t)) in distribution on Z)([0, 1]). 
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4. Moderate deviations for GEM. Let a{9) satisfy (3.5). 

The MDP for GEM is thus the LDP for the family of {a{e){Xi{e),X2{e), . . .) : 
6 > 0} when 6 approaches infinity. The result is proved through explicit cal- 
culations. 

Theorem 4.1. The family {P{a{9){Xi{e),X2{9), . . .) e ■) -.9 > 0} satis- 
fies a LDP on with speed and rate function 

oo 

(4.1) I(xi,X2, . . .) = ^Xi. 

i=l 

Proof. Let us first prove the local LDP. For any x,y in , set 

\xi - yil A 1 



|x-y| =^ 



2' 



For any x in MP^ and any 6 > 0, one can choose n sufficiently large that 

2» 



Y.Zn h < ^1'^- ^^^'^ ^1 < we have 



(4.2) (yGM^: max - x,| < ^ij C {y G : |y - x| < 5}. 

l<i<n J 

By taking limits in (4.2), in the order 9 oo,(5i — > 0,n ^ cx),(5 — > 0, it 
follows that 

limliminf 4?^logP(|a(0)(Xi(0),X2(0),...) -x| <5) 

(4.3) 

> lim lim liminf logpf max |a(0)Xj(0) — Xjl < 5i 
On the other hand, for any m > 1 such that 2"*" > b 

{yGR^:|y-x|<<^}c(yGM^: max |yi - x^j < 2™5l. 

l<i<m J 

Thus for any (52 < 5, 

(4.4) {yGM^:|y-x|<<52}c(yGM?': max \yi- x,\<2'^b\. 

y l<i<m ) 

By taking the limits in (4.4), in the order 9 — > oo, 62 ^ 0, 6 ^ 0, m ^ 00, 
it follows that 

lim limsup^logP(|a(0)(Xi(^),X2(e),...)-x| <52) 

(4.5) < lim lim lim sup ^^^^ 

X logPf max |a((9)Xi(6') - Xi\ < 2™(5 ) . 

\l<i<m / 
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It is known (cf. page 107 of [1]) that for J27=i Hi where < < l,k 
1, . . . ,n, the joint density function of {Xi{9), . . . ,Xn{9)) is 



(1 - yi)(l - (yi + y2)){l - (yi + • • • + Vn-i)) ' 
For any n > 1, (5 > 0, it follows from (4.6) that for sufficiently large 9 

xi H hx„ + n(5\^-V SO 



X 1 



a{e) J \a{9) 

— S\ f xi + X2 — 2(5 



X 1 



a{e) J\ a{9) 
xi H l-Xn-i - (n - 1)5 



< P max \a{e)Xi{e) -x^\<6 

\ l<i<n 

Kill - ~ n5\^~^ 



a{9) J \a{9) 

, , xi + 6\( xi + X2 + 2(5\ 

Xl^ h Xn-i + {n- 1)5 

W) 



Therefore 



lim lim sup — ^ log P ( max \a{9)Xi{9) — Xi\ <5 

= lim lim inf ^ log P ( max |a(0)Xj(0) — Xi\ < 5 

n 
i=l 

which combined with (4.3) and (4.5) implies that 

lim lim sup ^ logP{\a{9){Xi{9),X2{9), . . .) - x| < 5) 

= lim liminf ^ logP{\a{9){Xi{9), X2{9), ...)-x.\<5) 

oo 

1=1 



MODERATE DEVIATIONS 



13 



Now we show the exponential tightness. For any n > 1 and L > 1, it follows 
from direct calculation that 

p{a{e)Xn{e) >L)< P{a{e)Un >l) = (i ^ 



a{e) 



where (1- ^)+ is the positive part of (1-^). Set K = YY^^[Q,iL\. Then 
X is a compact subset of . Noting that for a; > 

(l-x)+<e-^ 

we get 

oo 

p{a{e){Xi{e),X2{e), ...)iK) <Y,p{a{e)Xi{9) > n) 

i=l 

oo / • r 



a(9) 



i=l ^ 



which implies 



limsup^logP(a(^)(Xi(0),X2(e),...) ^K) < -L. 

5. Moderate deviations for the Poisson Dirichlet distribution. Theo- 
rem 3.3 says that F{0) = {Pi{9) , P2{9) , ■ ■ .) approaches a nontrivial random 
sequence when scaled by a factor of 6 and shifted by P{6). Replacing the 
scaling factor by a{6) satisfying (3.5), we get 

(5.1) a(0)(^p(0)_M(i,i,...))^(O,o,...). 

The LDP corresponds to the case when a{9) = 1 and has been established 
in [4]. In this section, we establish the MDP for P{e) = {Pi{e), P2{9), . . .) 
or, equivalently, the LDP associated with the limits in (5.1). Considering 
the connection to Poisson point process, it is thus natural to start with the 
MDP for i(ai(0),a2(e),...). 

We first establish the MDP for an{9)/9 for any n followed by the MDP 
for {ai{9)/9,. . . ,an{9)/9). The infinite-dimensional case follows from finite- 
dimensional approximation. To go from the MDP for ^{ai{9),a2{9), . . .) to 
the MDP for the Poisson-Dirichlet distribution, one would hope to prove 
that a certain exponential equivalency holds. 
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5.1. MDP for "^"g^^ . It is known (cf. [15]) that for each n > 1, the density 
function of {ai(9), . . . ,an{0)) is 

(5.2) fniui,...,Un)= ^" e-Sr=i"'-^^i(""), ni>n2>--->u,>0. 

In particular, the density function of cri{9) is 

(5.3) 0u-ie-"-^^i("), M>0, 

where Ei{u) = y~^e~y dy. We extend Ei{u) to the whole real line by 
defining Ei(u) = +oo for n < 0. 

The distribution function of o"i(^) is 

(5.4) P{(Ji(0)<n} = 6-^-^1 u>0. 

One can find on page 146 in [15] the following explicit expression for the 
distribution function of cr„(^) for all n > 1. 

Lemma 5.1. The distribution function of an{0) is 

1 

(5.5) Fn{y) = -. — / n"-ie-" du, y>0. 

Next we establish the MDP for ai{9)/e. 

Theorem 5.1. The MDP holds for ai{9)/e with speed ^ and rate 
function 



Ji{x) 



X, x>0, 
oo, otherwise. 



Proof. For any fixed x, we have 

(5.6) p|a(0)(^^li^l^:M^ =e-ei^i(W«W)-+/3W). 

By L'Hospital's rule, 

(5.7) lim xe^'Eiix) = lim = 1. 

Restricting to a subsequence if necessary we can assume without loss of 
generality that lim^^oo ['^^^ + Pi^^))] exists in [— oo, +oo] . If the limit is neg- 
ative, then the event {a(0)(^^-^^t^^) < x} is eventually empty. Therefore 



(5.8) li^_^sup^logp(aw(""'>/"")<. 



-OO. 
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If lim0_»oo[^^2; + [i{9)] is a nonnegative finite number, then x is negative 
and -^Q^ and [5{9) are of the same scale as \og9. It follows from (5.6) that 
(5.8) also holds in this case. 

When lime^oo[^^3; + /3(^)] = oo, we can use (5.7) to get 



lim 

0-»oo 



a{e) 



(5.9) 



lim 



a{9) 



a{9) 



x + (3{9) 
log 9 



-{e/a(8))x 



oo 9 [9/a{9))x + (3{9) 



0, 



-oo, 



Thus 



(5.10) lim ^logP(a(0) 



X > 
X < 0. 



oi{9)-m 



<x =0, 



and 



a{9) 



(5.11) limsup^logP 0(6*) 



a^{9)-m \ 

9 y 



< X 



x>0, 



-00, X < 0. 



For X > 0, it follows from (5.6) and (5.7) that 

■a^{9)-m 



(5.12) 



lim sup log P[a{9) 



limsup — — log[l — e 

e-*oc 9 



> X 



r «(^) 1 
iim sup — — log 

0^00 9 



9E^ 



eEi{{a{e)/e)x+[3{e)) 
a{9) 



-x + (5{9) 



< —X. 



Together, (5.12) and (5.11) imply that the family of the laws of a(^)(^^-^^^g-^^~ 
is exponentially tight. 

Let gi{u) denote the density function of a{9){^-^^^^^^^^) . Then it follows 
from (5.3) that 

6» log6' 
(5.13) + 

This, combined with (5.9), implies that 
a{9) 



^^{e/a{e))u^~eE^{{e/a{e))u+m)_ 



(5.14) 
(5.15) 



9 
a{9) 



\oggi{u)^-u, 'u>0, 
log5i(ii) — > — 00, < 0. 
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For each 2; 7^ 0, choose 6 smah enough so that all numbers in the interval 
[x — 5, X + (5] are of the same sign. It is not hard to see that for u£ [x — 5,x + 6], 

loge 



9i{u) > 



(5.16) 
and 



a{e) {e/a{e)){x + 5) + p{e) 

^ ^^{e/a(e))(x+5)^~eE^({e/a{e))(x-&)+m) 



e \oge 

9i[u) < 



(5.17) 



aiO) {e/a{e)){x-5)+l3{e) 
^ ^-{e/a{e)){x-S)^-eEr{{e/a{e)){x+S)+m) _ 
Putting (5.14), (5.15), (5.16) and (5.17) together, we get that for x > 0, 
hm Mm sup ^ log P (a{e) ( ^ll^l^M) g _ 5, ^ + 5)) 

= lim liminf ^ log pfa(0) ( ^M^M.) ^ _ 5, ^ + 5) 



X, 



lim lim sup log P f a(^) f ^i^^V^) E (rr - X + (5) ) = -00. 



and for any x < 

Together, (5.10) and (5.11) imply that 

lim lim sup ^ log P(a{e) ( ^^^^^"^^^^ ) G (-5, 6) 



>oo 



lim liminf 4^ logpfa(^)f^l^^l^M U ) =0. 



The theorem now follows from the local LDP and exponential tightness. 

□ 

The next theorem gives the MDP of an{6)/6 for n > 2. 

Theorem 5.2. The MDP holds for an{0)/O with speed ^ and rate 
function Jn{x) = nx, x > 0. 

Proof. For x > 0, it follows from (5.10) that 

„,,^^^.,P(,,(«)),.) 
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By L'Hospital's rule, 



(5.19) limny-"e^ ru"-^e-"(in = l. 

y^o Jo 

Thus it follows from Lemma 5.1, (5.9) and (5.19) that 

..^,„,P(„<,(.M)),.) 

9^oo 6 

+ lim ^nlogf eEi( -^x + (3(9) 



(5.20) 

= n lim ^logf e-(^/«W)- 

""eToo ^[{9/ai9))x + p{0) 

= —nx. 

For X < 0, as in the proof of Theorem 5.1, it suffices to obtain estimates 
for those x such that 



lim -j^x + P{0) =+oo. 

Since 0Ei{^x + p{9)) « (^e/am'+m ^'^'^''^'^^'' approaches infinity as 
tends to infinity, one gets that 

limsup^logp(a(^)(-"(^)/(^))<x 

a{0) 

P 

(5.21) 



limsup ^ log ^e-'^ima{e)).+m) 



-oo. 



The exponential tightness of the laws of {a{0) """(^)~^(^) | j^qw follows from 
(5.20) and (5.21). The local LDP can be obtained by an argument similar 
to that used in Theorem 5.1. □ 

5.2. MDPfor Uai{0),a2{0), ■ . .). For each 7i > 2, we have 
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Theorem 5.3. The family {P{a{e) C'%'^^'^'^ '""W-^(^) ) e-):9> 
0} satisfies a LDP on M" with speed and rate function 



(5.22) /„( 



^Xi, ifO<Xn<---<Xi, 

i=l 

+00, otherwise. 



Proof. It follows from (5.2) that for xi > X2 - ■ ■ > Xn and -^^Xn + 
P{e) > 0, the density function of ^{(7ii9) - f3{e), . . . ,ct„(0) - ^{6)) is 

log(^) \ 



/ \ " / " 



i6))xi+m 

(5.23) 

X ^-[ie/ai9))J2lU^^+0Eli(e/a(e))x„+m)]_ 

By direct calculation, 

a(6) " 

(5.24) —^loggn{xi,...,Xn)^-'^Xi, Xn>0, 

1=1 

(5.25) ^^^loggn{xi,...,Xn) ^ -oo, Xn<0. 

For xi > X2 • • • > let -B((xi, . . . ,Xn),S) denote the closed ball centered 
at {xi, . . . ,Xn) with radius 5, and B°{{xi, . . . ,Xn),S) be the corresponding 
open ball. Then for x.„ > 0, 

hmlimsup ^ logpf ^(ai(e) - m, ...,M0)- m) 

£ B{{xi,. . . ,Xn),5) 

(5.26) = hm hminf ^ fogpf 4r^(ai(0) - /3(^), . . . , a„(e) - P{9)) 



5->o e- 



GB°((xi,...,2;J,<5) 



= -E 

and for any x„ < 



^2 

j = l 



lim limsup ^ log p(^{ai{9) - P{9), . . . , a„(0) - /3(0)) 



5^0 



eB{{xi,...,Xn),S)^ 
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(5.27) =Mmmf^logp(^{ai{e)-m,---,MO)-m) 



5^0 e^oo 



-oo. 



If xi = 0, the upper estimate follows from Theorem 5.1. If Xr-i > 0,Xr = 
for some 1 < r < n, then the upper estimate is obtained from that of 
^^((Ti(0) — P{9),. . . ,ar-i{9) — P{9)). The lower estimate when Xr = for 
some 1 < r < n is obtained by approximating the boundary with open sub- 
sets that have all positive coordinates. 

Fix an L > 0. Noting that [j7=i{^{M0) - l3{0)) > L} = {^(cJi(^) - 
P{9)) > L}, it follows that 

(5.28) hm hmsup^logp|u(^(^..(^)-/3(^))>^|| = -oo. 



ia,{9)-m)<-L 



On the other hand, 

hmsup— -logP<^ Ul 
(5.29) 

< limsup^logp|^(a„(0) - m) < -^} = -oo. 
Therefore we have the exponential tightness and the theorem. □ 
The MDP for ^{ai{9),a2{9), . . .) is derived in the next theorem. 

Theorem 5.4. The family {P(^ (di (9) - p{9),a2{9) - /3{9), ...)£■): 9 > 
0} satisfies a LDP on with speed and rate function 



(5.30) I{xi,X2, 



oo 

^Xi, xi>--->0, 

i=l 

OO, otherwise. 



Proof. Identify M°° with the projective limit of M",n = 1, Then 

the theorem follows from Theorem 3.3 in [5] and Theorem 5.3. □ 

5.3. MDP for the Poisson-Dirichlet distribution. Using the results in 
the previous subsection we now derive the MDP for the Poisson-Dirichlet 
distribution. The representation (1.1), combined with the fact that cj(0) is 
approximately 9, seems to suggest that the MDP for the Poisson-Dirichlet 
distribution should follow from the MDP for h{(Ji{0),a2{9), . . .). This turns 
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out to be true. It does not seem to be easy to get a more direct proof using 
the explicit expression in [28] of the density functions of {Pi{9), . . . ,Pn{6)) 
for each n > 1. 

Theorem 5.5. For each n > 1, the family {P{a{9){Pi{9) - 

Pn{0) - ^,...) G-):e>0} satisfies a LDP on ]R°° with speed ^ and 
rate function 



(5.31) 



I{xi,X2, 



^Xi, xi>--->0, 
otherwise. 



i=l 

oo. 



Proof. From representation (1.1), one obtains that 



(5.32) 



aiO){PniO)-^) 



a{e) 



a{e) 



cTn{e)-m 



Write 



a{e)m 



and without loss of generality we assume that 



lim 7(0) = cG [0,+oo]. 

S— >oo 



It is clear that 



(5.33) 



a{e) 



oo, — > oo. 



If c < oo, it follows from Corollary 3.1 that for any L > 

a{e) 



(5.34) limsup^logp(7(6l) 
For c = oo, and any 1 > (5 > 



a{e) 



1 



> L} = -oo. 
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m 

rem 3.2) for a{e)/e, and (5.33) that 



Since 'y{9) < fj{6) and lim^^ 



: 0, it fohows from the MDP (Theo- 



(5.36) 



hmsup^^^ log-P<{ 7(0) 



a{e) 



>{l-5)L^ 



> (1-<5)L 



-00, 



which, combined with Corollary 3.1 and (5.35), shows that (5.34) still holds 
in this case. Therefore a{e){Pn{e) - ^) and ^a{9)[^^^^^^^] are expo- 
nentially equivalent. 

Since ^a(g)[ '^"<%~^^^) ] is exponentially equivalent to a(g)[ '^"^%~^^^) ] 

by Lemma 2.1 and Corollary 3.1, it follows that a{9){Pn{9) — ^^) and 
^|-g^|- o-n(6')-/3(e) j g^^g exponentially equivalent for all n > 1. Thus the MDP for 
the Poisson-Dirichlet distribution is the same as the MDP for ^((Ti(0),o"2(0), 
...). □ 



6. Moderate deviations for homozygosity. For each m > 2, it is shown 
in Theorem 3.4 that the scaled homozygosity p^^^ Hm(P{d)) satisfies a 

fluctuation theorem with c{6) = 1 and b{9) = VO. It is thus natural to 
consider MDPs for p"^^ Hm(P{0)) or equivalently the LDP for the family 

{a{e)[^^Hm(P{9)) - 1] :0 > 0} for a scale a{e) satisfying 

a{9) 



(6.1) 



lim a{9) = 00, 



lim ^ 



0. 



It will turn out in Remark 6.1 that the following additional restriction on 
a{9) is necessary in order to get the MDP: for some < e < l/(2m — 1), 



(6.2) 



,1-e 



liminf —, — 

e~*oo 6)(m.~l)/{2m.-l) 



>0. 



The main idea of the proof is to explore the connection between homozy- 
gosity and the Poisson process, and apply Campbell's theorem. 
Let us first consider the MDP of 



G 



-Y.a]\9). 



It follows from Campbell's theorem that 



'dy 



22 



S. FENG AND F. GAO 



which imphes that {G'f\e> 0} is a random process with stationary and 
independent increments. The difficulty here is that the exponential moment 
is not finite. MDPs for models with infinite exponential moment have been 
studied in [14, 17, 22]. A typical way of establishing the MDP is to verify 
the following Ledoux condition [22] : there exists a constant M > such that 
for any 6 > 0, 



This condition does not seem to be easy to verify for our model. Therefore 
we employ a truncation procedure. 

Lemma 6.1. Set 

Ge = {a{e)-e,G'^"'^ -T{m)e). 

Then the family {^^^Gq-.O > 0} satisfies a LDP with speed " and rate 
function 

A*(x,y) := \ {V{2m)x^ - 2V{m + l)xy + y^), 

2(1 (2mj — i (m + \y) 

Proof. By (6.1) and (6.2), there exist r > and a positive integer 
i > 3 V (2m-i)£ such that 

lim — — = +00 

and 



Ql-2/({2m-l)0(6)) 

lim -, -YTTT^ n — = OO- 

e-,oo gi'm-l)/{2m.-l) 



Take 

Q(/-2)/{m-l)Z(^) 



7(0) 



log((a2(^)/^)aa-2)/(m-i)«(0))- 
Then '~f{0) grows faster than a positive power of 9 and 

-f(9) , a'^idhid) 

lim ,, „, TTTTT = 0, lim = oo. 

e^oo a(^-2)/(™-i)^(6') e^oo 9 



MODERATE DEVIATIONS 



23 



Set 



and 



Define 

A(a, /3) = \{a^ + 2T{m + l)a/3 + T{2m)l3'^) 

1 r(m + 1)\ /a^ 



r(m + i) r(2m) y V/3 



Then 



sup {ax + Py — A(a, (3)} 



(x y) 



r(2m) -r(m+ 1)\ (x 



2(r(2m)-r(m + l)2) ^ ^ ' \-T{m+l) 1 
-{T{2m)x^ - 2T{m + l)xy + y"^ 



2(r(2m) -r(m + i)2; 



For any a G R, /? G M, 



1 



logEl exp^ -77t(«(G^'^ - ^(G?^)) + P{G^r^ - E{Gt^))) 



a{e) 



A(a,/3) 



log exp|^ j 



,(l/a(e))(aj/+/3y'") 



■7(e) 



a2(^) / (e(iMe)){"s/+/33/'-) _ ^ _ a-i(0)(ay + /?y™))y"^e"^ dy 



A(«,/3) 



< 



7(e) 1 



{ay + l3y'^fy-^e-ydy-K{a,l3) 
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oo 1 

+ E T-,a'^''^\e){\a\ + \[5\^{er-'fm 

k=l+l 



as 6* — > oo. 



Therefore, by the Gartner-Elhs theorem, {^^^Gq, " , A*) satisfies a LDP. 
Noting that ^{6) grows faster than 0" for a certain a > 0, it follows that 

lim eEi{j{e)) = lim 07(0)-ie-T(^) =0. 

9— >oo B^oo 

Taking into account the fact that G^p < cr{9), G^^ < G^p\ we have that 
hm sup 4^ log P f I (a(^) , ) - (G« , ) | > 5 



< lim sup ^ log P (^1 ((7(0) , Gp^ ) |/|,, (e)>,(e)} > 



< lim sup log P(ai (0) > 7(0)) 

= liinsup^log(l-e-^^iW^))) 
^ J [l°g'' + logEl(T(9))l°"W 



= — CXD 



which implies that for any 5 > 0, 

u.n.np ^g) logpf IQ" - g° - (^(6^") i^(Qr') - r(>n)^)l , , 

(6.3) 

= — oo. 
By direct calculation, 

{E{GP),E{Gt^))-{e,T{ra)e) 
e/a{9) 

= - lim a{9)( / e^ydy, / /'""^e-^'dy 
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= - lim a(0)(e-^(^),72'"-i(e)e^^(^)) = O, 
6^00 

which, combined with (6.3), imphes that ^^^Gg and ^^^Gq are exponentially 
equivalent. Therefore 

satisfies the LDP. □ 

Now we are ready to prove the main result of this section. 

Theorem 6.1. The family a{9){Y^Hm{^{9))-l) satisfies a LDP with 



speed " i^^ and rate function 



2 



2(r{2m)/r(m)^-m^) ' 



Proof. By direct calculation, 
\L [m) 

0m— 1 



amu^] -1)+' * ^"'owiGr-Eior^ 



(^{0); J \a{e)j T{m)e 



"^'^?A^^ ^V.W; Tim) 



Noting that for any i>l and for any 6 > 0, 



lim log P 

e->oo 6 



1 



>5 =-oo. 



It then follows that 

~,(m) ,„(m) 

/ 1/ \ KII/PII nil 

are exponentially equivalent to 

a{e)m{e-a{e)) a{0){Gt'^ -E{G^^^)) 

respectively. Thus 

3771—1 



\i (mj 
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and 

a{e)m{e-a{0)) a{e){G^^^ - E{G^^'^)) 
9 T{m)9 

have the same LDP. 
Since 

iy/r{rny)-mx=z^ = 2(r(2m)/r(m)2 - m2) ' 

Lemma 6.1 and the contraction principle yield that 

a{e)m{e - a{9)) ^ a{e){G^^^ - £{0^"^)) a^O) 



T{m)e ' e ' 2(r(2m)/r(m)2 -m2) 

satisfies a LDP, and the theorem follows. □ 

Remark 6.1. Choose the scaling factor, a{9) = 6^ . Then the MDP ob- 
tained here requires that 7 lies between 2ni-i ^'^'^ \' natural to ask 
what happens for 7 < 2rn-i • follows from Lemma 6.1 and the contraction 

principle that the family { "^^^^'^'^ — — ^ : > 0} satisfies a LDP with 
speed and a rate function J{x) = 2r(2m) • Thus for any 5 > 0, there 

exists ^0 > such that for all 9 >8q, 

p(ior'-£(Gr)i>^)<e.p{-^[/TO-i/2i}. 

bmce is a random process with stationary and independent 

increments, one can find sufficiently small 61,62 > such that 

p(,cS™)-.(cS-)),>|^) 

<P( |gK - > +P( KG^ - EiG^^))] > 



a(^ + l); ' V ' " a{9) 



< 2exp 



A52)-l/2]|. 



o2(0)' 

The fact that lim2,.^oo Jix) = +00, yields 

limsuplimsup^logpf {G^^^ - E{Gt^)\ > ^) = -00 
5^00 0^00 9 \ a{9)J 

which, combined with the fact that E{G^^) is a finite number, implies 

(6.4) limsuplimsup — ^— ^ log pIg^^^ > ^ = —00. 

5^00 0^00 9 V a{9)J 
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Since 

it follows from (6.4) that 

limsuplimsup — - — log£/i — — = — oo. 

s^oo e^oo \\a[6)J J 

Using the relation (5.7), one gets 

m — 1 



7> 



2m- 1' 



which corresponds to the critical case of e = in (6.2). Thus the range of 
scaling obtained here is the best that one can get for the MDP with speed 

e ■ 

7. Moderate deviations for the Dirichlet process. In this section, the 
MDP for the Dirichlet process is derived through a combination of the LDP 
for the gamma distribution and MDPs for processes with stationary inde- 
pendent increments. 

The Dirichlet (9, z/) distribution can be represented by 

X{9u{[0,t])) 
X{6) ' 

where {X{t),t S [0,oo)} is a Gamma process. By Theorem 3.1, the family 
{P{X{0)/9 e •), ^ > 0} satisfies a LDP in R+ with speed 9 and rate function 
I{x) given by (3.4). 

Let a{9) be a positive function satisfying (6.1). With a time deformation, 
the following theorem is a minor generalization of the result in [25]. For 
completeness, a sketched proof is included. 

Theorem 7.1. Let {^(t),t G [0,oo)} he a stochastic process with sta- 
tionary independent increments and right- continuous paths with ^(0) = 0, 
^(^(1)) = 1, Var(e(l)) = I, and 

£;(e5|C{i)l ) < oo, for some 6 > 0. 
Let u he a finite measure on [0, 1] such that ^{{t}) = for all t G [0, 1] . Define 

and 

ieAt)=ieAM). v{t) = v{[Q,t\). 
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Then the family {P{a{e){^eAt) " ^^(0)/^ e-):9>0} satisfies a LDP in 
(D[0, 1], II • II) with speed and rate function 

I +CXD, otherwise, 
where \\ip\\ := sup^ejo,!] |93(i)l for ip € D[0,1]. 

Proof. It suffices to verify the following three conclusions (cf. [2, 31]): 
(i) For any < ti < t2 < • • • < ifc < 1, 

satisfies a LDP with the speed and the rate function 



1 ^ I _ |2 

(ii) For any 5 > 0, 

lim sup limsup — - — logPl sup |^(0z/((s, t])) 



e^Oj,g[o,l] e^oo d \s<t<s+e 

(7.1) 

-M(M])l>^^)=-oo. 
(iii) I{ip) = suptj^...^j^c(o,i] hu-.tdAh), Atk))- 

Since ^(t) is a random process with stationary and independent incre- 
ments and the mapping: 

{ZI,Z2- Zl,...,Zk- Zk-l) {zi,Z2,...,Zk) 

is continuous in M^, it is easy to get (i) from the Gartner-Ellis theorem, the 
product principle and the contraction principle; (iii) is a consequence of the 
Cauchy-Schwarz inequality and the submartingale convergence theorem. 

Finally, we verify (ii). By Corollary 4 in [24], it is easy to see that there 
is a universal constant c > 1 such that 

p( sup ie(^z.((s,t]))-M(.,t])i> ^ 

65 



s<t<s+e a{6) J 

<cp(\i{ev{{s,s + e]))-ev{{s,s + e])\ > 



a{e)c 



exp|-i^a|e(0z/((5, s + e])) - eu{{s, s + e])||) 
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a(e(i) - 



-1 



a{9) 



a(e(l)-l) 



eu{{s,s+e]) 



eu({s,s+e]) 



where a > is arbitrary. By the hypotheses, expanding the cumulant yields 

r a2 



2a2(0) 



o 1 



Therefore 



hm sup 



a\e) 



\ogP[ sup \^{eu{{s,t\))-Oy{{sM>^.^ 

s<t<s+e CL[U) 



( a6 a^uUsjS + e]) 
< -sup\ 



2c^iy{{s,s + e]) 



which imphes (7.1). □ 

We now estabUsh the MDP for the Dirichlet process. 

Theorem 7.2. The family {P{a{e){Z0^,y{t) - v{t)) G •) > 0} satisfies 
a LDP in D[0, 1] with speed —4-^ and rate function 



'D(f) 



2 Jo 

+ 00, 



dip 
dv 



it) 



u{dt), ifLp<^v, =0, 
otherwise. 



Proof. Choose ^(i) in Theorem 7.1 to be the Gamma process X{t). 



Set 



Ye At) 



a{e){XeAt) - Oi^it)) a{e)y{t){XeA'^) - 9) 



(l-uit)) 



a{9){XeAt) - OHt)) 



u{t) 



aie)iXeA'^)-XeAt)-0{l-iym 



By Theorem 7.1 and the contraction principle, the family {P(Yg^^{t) £ ■) -.O > 
0} satisfies a LDP in Z)[0, 1] with speed " and rate function 

inf{/(V); - A^)Ht) = V{t),t G [0, 1]} 



— inf 

2 aeK 
+ 00, 



i^{dt), if (^(1) = 0, 

otherwise, 



30 



S. FENG AND F. GAO 



1 

2 Jo 

+ 00, 



dip 
dv 



u{dt), if 99 < I/, (^(1) = 0, 
otherwise. 



Since 

\a{e){Ze,u{t)-u{t))-Ye,,{t)\ 



e 



a{e){Xe,u{t)-0^{t)) 



it follows that for any (5 > 0, and e > 0, 

p( sup \a{e){Ze,,{t) - u{t)) - Ye^,{t)\ > ^ 

Vi6fO,ll 



+ 



a{e)v{t){XeA'^)-e) 



l[0 

< p 



1 



+ p 



a{9){Xe,,a) - 6) 



> e ) + P[ sup 

ie[o,i] 



ai9)iXe,u{t) - OHt)) 



> 



2e 



> 



2e 



Now from the LDP of Xe^^(l), one obtains 
1 



lim sup - log P 
which implies 



1 



>e]<- inf I{x) < 0, 

' |l/2:-l|>£ 



lim sup — - — log P 
From the MDP of Xq. we have 

lim lim sup — - — log P I sup 

^ Vte[o,i] 



Xa,A\^ 



> e 



a{e){XeAt)-Gy{t)) 



-00. 



and 



lim lim sup — -- — log P 
Therefore, for any 5 > 0, 
lim sup — - — 



a{e){XeA'^)-e) 



>-|=-oo. 



logP( sup \a{e){ZeAt) - ^^{t)) - > 5 

te[o,i] 



-00; 



that is, a{9){ZQy{t) — fit)) is exponentially equivalent to 10,1/(0 • ^ 
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